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We present the Lagrangian and Hamiltonian framework which incorporates null dust as a source 
into canonical gravity. Null dust is a generalized Lagrangian system which is described by six 
Clebsch potentials of its four-velocity Pfaff form. The Dirac-ADM decomposition splits these into 
three canonical coordinates (the comoving coordinates of the dust) and their conjugate momenta 
(appropriate projections of four-velocity). Unlike ordinary dust of massive particles, null dust 
' therefore has three rather than four degrees of freedom per space point. These are evolved by a 

Hamiltonian which is a linear combination of energy and momentum densities of the dust. The 
energy density is the norm of the momentum density with respect to the spatial metric. The 
' coupling to geometry is achieved by adding these densities to the gravitational super-Hamiltonian 

.^fH , and supermomentum. This leads to appropriate Hamiltonian and momentum constraints in the 

■ phase space of the system. The constraints can be rewritten in two alternative forms in which they 

00 ' generate a true Lie algebra. The Dirac constraint quantization of the system is formally accomplished 

by imposing the new constraints as quantum operator restrictions on state functionals. We compare 
the canonical schemes for null and ordinary dust and emhasize their differences. 
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' Null dust has been widely used as a simple matter source both in classical and semiclassical gravity. Its equations 
^ of motion follow directly from the conservation of the energy-momentum tensor. However, the inclusion of null dust 
as a source into canonical gravity requires careful identification of its own dynamical degrees of freedom. For this 
I ■ purpose, one needs to construct a spacetime action depending on appropriate Eulerian variables and bring it into 
bi). canonical form by the Dirac-ADM (Arnowitt, Deser and Misner) procedure. The coupling to gravity, like that of 
other nonderivative systems, is then entirely straightforward. The ordinary dust of massive particles was treated in 
this manner by Brown and Kuchai^ 0|. Our goal is to develop a similar formalism for null dust. 
' The main application which we have in mind is minisuperspace and midisuperspace quantization of canonical models 
^ , which include null dust as a source. The specific models based on null dust are both numerous and simple. After the 
- - ' discovery of Vaidya's 'radiating Schwarzschild metric' there were found many other, more general exact solutions 
of Einstein's equations with null dust as a matter source. Above all, such models have recently been used to clarify the 
formation of naked singularities during a spherical gravitational collapse, to describe mass inflation inside black holes, 
and to model the formation and Hawking evaporation of black holes. We briefly review these topics in Appendix B. 
Our formalism is designed for studying such issues in quantum rather than in classical or semiclassical contexts. 

Null dust is intimately connected with the behavior of zero-rest-mass fields in geometrical optics limit. The 
energy-momentum tensor of such fields takes in that limit the form of the energy- momentum tensor of null dust. 
One can then reinterpret some exact solutions of Einstein's equations with null dust as spacetimes produced by 
zero-rest-mass (in particular, electromagnetic) fields. Careful studies of the high-frequency limit of the gravitational 
radiation itself revealed that it also can be described by the energy-momentum tensor of null dust. Moreover, which 
is especially relevant for the present paper, such a connection can be established at the level of a variational principle. 
All of this indicates that null dust is much more closely related to fundamental fields than ordinary dust formed by 
phenomenological massive particles. We explain some of these connections in Appendix A. 



I. INTRODUCTION 
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We start our exposition by reviewing how the dynamics of incoherent dust follows from the conservation law of the 
energy-momentum tensor (Section 2). This enables us to pinpoint at the very beginning the main difference between 
ordinary dust and null dust: The normalization of timelike four-velocity selects its parametrization by proper time, 
the null normalization of lightlikc four-velocity leaves its parametrization arbitrary. This is of paramount importance 
both for the Lagrangian and Hamiltonian descriptions of null dust. 

Since null geodesies are somewhat less intuitive than timelike geodesies, we briefly summarize the basic properties 
of null congruences in Section 3. We explain how to obtain an affine parametrization of such congruences, but stress 
its essential ambiguity which prevents the unique separation of mass distribution of null dust particles from their 
four-velocity. When it comes to producing the gravitational field, the mass distribution can simply be reabsorbed 
into four-velocity, which is the reason why it does not naturally occur as a separate variable either in the Lagrangian 
or Hamiltonian frameworks. We formulate a spacetime variational principle from which the null geodesic equations 
of motion follow in Section 4. The variational principle for null dust is quite similar to the variational principle for 
ordinary dust given by Brown and Kuchaf but there are several characteristic differences. The most important 
one was already mentioned: Null worldlines have no natural parametrization and hence the null velocity appears 
in the variational principle as the Pfaff form of six scalars (the Clebsch potentials) rather than seven scalars which 
characterize the timelike velocity. Consequences of this distinction can be traced throughout the whole formalism. 
In Appendix C we illustrate on explicit examples the decomposition of null covector fields into Clebsch potentials 
which is a prerequisite of our variational principle. In Section 5 we show that any solution of the Euler equations 
of the variational principle provides enough building blocks to reconstruct an afhnely parametrized four-velocity fc" 
and the associated mass distribution M. We also discuss other special parametrizations. In Section 6, we cast our 
covariant spacetime action into Hamiltonian form by following the Dirac-ADM algorithm. The details of this process 
are substantially different from the steps which need to be taken for ordinary dust. We express the energy and 
momentum densities of null dust in terms of appropriate canonical variables. The energy density turns out to be the 
norm of the momentum density with respect to the spatial metric. It transpires that null dust has only three degrees 
of freedom per space point, one less than ordinary dust which has four. The missing degree of freedom is a privileged 
scalar parameter (like proper time) along lightlike geodesies. The missing canonically conjugate momentum is the 
mass distribution which has been reabsorbed into the four-velocity form. We conclude this section by writing down 
the standard Hamiltonian and momentum constraints for geometry coupled to null dust. In Section 7, we rewrite 
these constraints in two alternative forms in which they generate a true Lie algebra. In this process, the Hamiltonian 
constraint is replaced by alternative constraints which contain only geometric variables. This feature of the constraint 
is related to a Rainich-type 'already unified theory' for geometry coupled to null dust. In Section 8, we show how one 
can formally impose the new form of constraints as quantum operator restrictions on state functionals. The outcome 
of this procedure is a single functional differential equation for physical state functionals ^ [g] which depend solely on 
the spatial metric g in the dust frame. In the final Section 9, we compare the canonical formalism and the ensuing 
quantum theory for null dust with those for ordinary dust and emphasize their differences. 

Our conventions follow those of Misner, Thorne and Wheeler ||^ , except for our choice of units which are such that 
167rG = 1 = c. 



II. DUST AS A SOURCE OF GRAVITY 



Incoherent dust is one of the simplest phenomenological sources of gravity in general relativity. Its energy- 
momentum tensor 

T"^ = MU^U^, (2.1) 



U''Ua = -l, (2.2) 
curves the spacetime according to Einstein's law of gravitation 

Dust is described by the four-velocity C/" of its particles and the (rest) mass density M of their distribution. The 
equations of motion of the dust are entirely contained in the energy-momentum conservation law 

V/jT"'^ = (2.4) 
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which follows from the Einstein law (2.3) through the Bianchi identities. The structure (2.1) of the energy-momentum 
tensor allows us to write Eq. (2.4) in the form 

MU'^VffU"' +V0{MU'^)U" = 0. (2.5) 

One sees that 

MUf^VfjU" oc U", (2.6) 

i.e., that the dust particles move along geodesies. The normalization (2.2) of the four-velocity tells us that the particle 
worldlines are parametrized by proper time. When one multiplies the geodesic equation by Ua , the normalization 
condition (2.2) implies the rest mass conservation 

V0{MU'^)=O. (2.7) 

By using Eq. (2.7) back in Eq. (2.5) one learns that proper time is an afBne parameter: 

Ul^WfjU" = . (2.8) 

These facts describe in full detail the motion of ordinary dust of massive particles. 

Null dust has the same energy-momentum tensor (2.1) as ordinary dust, but its particles are assumed to follow 
lightlike worldlines: 

U°'Uc=0. (2.9) 

The energy-momentum conservation law (2.4) still implies that those worldlines are geodesies, Eqs. (2.5)-(2.6). 
However, the null normalization (2.9) no longer enforces either the conservation law (2.7) or afRne parametrization 
(2.8). 

For ordinary dust, the decomposition of the energy-momentum tensor into the mass density M and four- velocity 
Ua is unique due to the timelike normalization (2.2). For null dust, the lightlike normalization (2.9) is preserved by 
an arbitrary scaling of Ua ■ 

U°'=AU°', A>0. (2.10) 

(The limitation A > is needed to preserve the future-pointing orientation of the worldlines.) By simultaneously 
rescaling the scalar M , 

M = A-2m, (2.11) 

one preserves the form (2.1) of the energy-momentum tensor. This shows that the decomposition of T"^ into M and 
U" is arbitrary. In particular, by taking A = M , one can eliminate the scalar M altogether and write T"^ entirely 
in terms of a single null vector 

/" := M^/'^U" (2.12) 

as 

rpaH ^ ;a;/3_ ^2.13) 

In terms of l", the geodesic equation (2.5) takes the form 

if^Vpl" + = 0. (2.14) 

We shall see later that this choice maximally simplifies the form of the null dust action and its canonical decomposition. 



III. NULL GEODESIC CONGRUENCES 



In a region of spacetime, A4, which is filled by dust whose worldlines do not intersect, the vector field U" defines 
a line congruence S. This congruence can be viewed as an abstract three-dimensional space, the 'dust space', whose 
points are the individual worldlines. The worldlines z G S can be locally labeled by three parameters z''{z) which 
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introduce a coordinate chart in S . We shall use the indices i, j, k from the middle of the Latin alphabet to denote 
the components of the objects in S ; they take the values 1,2,3. (A global standpoint replacing this local description 
is discussed in [0.) 

Through each event of the region there passes one and only one worldline. One can uniquely assign to each event 
y the labels z'^ of that worldline: 

z'' = Z'^{y). (3.1) 

Our interpretation of the scalar fields Z^{y) presupposes that their values constitute a good chart in S . Therefore, 
the three gradients must be three linearly independent covectors:|] 

U'- oc^S''^^'Z\pZ',-,Z\sS^jk^O. (3.2) 

Parametrize the curves of 5 by a parameter u whose rate of change judged by the size of U" is unity. In other words, 
if 

u = U{y) (3.3) 
is the value of u on the curve of S which passes through the event y, it holds that 

U°'V^U=l. (3.4) 
Equations (3.2) and (3.4) imply that Z^ — {U, Z^) are four independent functions of spacetime coordinates y": 

det(Z^,„) = i^"''^* U,^Z\pZ\^Z\s 5^Jk + . (3.5) 

The mapping Z : M. —* R S given locally by Eqs. (3.3) and (3.1) can thus be inverted into the mapping 
T : Ry. S ^ M. given locally by 

2/" = T"(^.,z'=). (3.6) 
Here, z^ distinguish different curves of the congruence S and u specifies the point on a given curve. The four vectors 

C/" = T",,, Z," = T",fc (3.7) 

form a basis in TM. dual to the cobasis 

^"^.a = ((7,c. , Z^„) (3.8) 

in T* M. ■ The basis (3.7) and cobasis (3.8) satisfy the standard orthonormality and completeness relations. In 
particular 

zfz^„-^t. (3.9) 

So far, everything applies equally well both to timelike and null congruences. Brown and Kuchaf ||l[ specialized 
the formalism to timelike congruences and applied it to Lagrangian description of ordinary dust. In this paper, we 
first briefly recapitulate how to specialize the formalism to null geodesic congruences (see, e.g., Q and H] for more 
details) and then use it for Lagrangian description of null dust. 

A geodesic null congruence J7" must satisfy the geodesic condition (2.6) and the null condition (2.9). These 
conditions still hold when the vector field U" is scaled by an arbitrary factor, Eq. (2.10). Instead of using that 
scaling for eliminating M from the energy-momentum tensor, Eqs. (2.11)-(2.13), one can use it for enforcing affine 
parametrization. In terms of the geodesic equation takes the form (2.14). Unless P happens to be divergencefree, 
it is not affinely parametrized. Let us first show that there exists a positive scaling factor 

A(y) = e^(^) (3.10) 



^The contravariant tensor density S"'^"'^ of weight 1 is the alternating symbol in M . The covariant tensor density Sijk of 
weight —1 is the alternating symbol in S . The Levi-Civita pseudotensor in Ai is denoted by e^'^T*. 
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such that 

V/3(A-i/^) =0. (3.11) 
The condition (3.11) amounts to a linear inhomogeneous equation 

l/^WpX = Vpl^ (3.12) 
for A . In the adapted coordinates u, z'^, Eq. (3.12) assumes the form 

'-^ = W^)iu,z). (3.13) 

Its general solution 

\{u,z)= [ du{W0l^){u,z) +Xoiz) (3.14) 
Jo 

depends on an arbitrary function Ao(2;) of z. By writing Eqs. (2.13) (2.14) and (3.11) in terms of the new variables 

k":=M°' and M:=K-^ (3.15) 
one learns that the vector field fc" is afHnely parametrized, 

fc/3V/3fc" = 0, (3.16) 
the mass distribution M satisfies the continuity equation 

V/3(MA;'') =0, (3.17) 
and the energy-momentum tensor takes the form 

T"^ = Mk'^k^. (3.18) 
The afEne parameter w is a monotonically increasing function of the old parameter u : 

v{u,z)= [ duA~^{u,z)+vo{z). (3.19) 



When we define a new mapping T^pp 

riMv,z)~r''{u{v,z),u) (3.20) 

we obtain 

fc« = ^^aff(^,^) _ (3 21) 

ov 

The afBne parametrization (3.19) depends on two arbitrary functions, Xo{z) and vo{z) . This means that along each 
geodesic the afHne parameter is determined only up to a linear transformation 

v = Ao^{z)v + vo, Ao{z)>0. (3.22) 

When we change the afJine parameter by Eq. (3.22), the null vector field fc"(y) is scaled into 

S"(y) = Ao(Z(y))fc"(2/). (3.23) 

The afHnely parametrized null vector field k°'{y) is thus determined only up to an arbitrary positive multiplicative 
factor Ao which is a function of comoving coordinates z'^ = Z''{y) . 

A congruence of affincly parametrized null geodesies is characterized by its twist (or rotation) w, expansion 9, and 
shear a. The corresponding scalars are given by 
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\ 1/2 

= I , (3.24) 

= \ (VaA:") , (3.25) 

-V^^kpjV'k^ -e^j , (3.26) 

where the square (round) brackets around indices denote antisymmetrization (symmetrization) . The twist can also 
be determined from the relation 

Luk'' = ie"^'^*fc^V5fc^ . (3.27) 

If ka is proportional to a gradient, 

fc«(y) = '/'(y)V',a(y), (3.28) 

the geodesies of S form null hypersurfaces ip — const to which fc" is orthogonal. A null geodesic congruence is 
hypersurface orthogonal if and only if it has a vanishing twist: uj — . 

Under the change (3.23) of affine parametrization, the rotation, expansion and shear all scale by the same factor: 

lo^Aqlo, e^Aae, Ao|cr| . (3.29) 

Also, by using Eqs. (3.10)-(3.12) and (3.15), one can reexpress them in terms of l" and its derivatives, and of the 
undifferentiated scaling factor (3.10), (3.14): 

/I 1 

u; = Af-V[„Z^]V"/'5 + -(V„r)M , (3.30) 

6i = A(V„r), (3.31) 

\a\ = A Qv(„Z^)V"?'5 - \{VJ^f^ ' . (3.32) 

The scalars (3.30)-(3.32) allow us to introduce other special parametrizations of null geodesic congruences. Consider- 
ations about the rate of expansion of a shadow image (see, e.g., Q and [^) lead to the concept of luminosity distance 
A. This is defined as any solution of the equation 

where v is an afhne parameter and 9 (assumed to be nonvanishing) is the expansion (3.25). Equations (3.15) and 
(3.23) ensure that the luminosity distance is identical with the scaling factor (3.10), (3.12). The luminosity distance 
played a prominent role in several classical works in radiation theory and in cosmology |^] . The mass distribution 
M introduced by Eq. (3.15) is the inverse square M = A^^ of the luminosity distance. The parallax distance p , 

p^e-\ (3.34) 

is also occasionally useful. 



IV. SPACETIME ACTION AND THE EULER EQUATIONS 

We describe null dust by six spacetime scalars Z*^, Wk ■ The interpretation of our state variables , Wk emerges 
from the form of the action and the resulting equations of motion. We shall see that Z'^ are comoving coordinates 
of null dust particles. By specifying the values z'" of the scalars Z^{y), we choose a particular null geodesic of the 
congruence S. The three gradients Z^ are assumed to be three independent covectors. We shall see later that none 
of them can be timelike. 

The four-velocity covector 1^ of a lightlike particle is given by its components Wk in the cobasis Z^ ■ 

L = WkZ\^. (4.1) 
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This relation expresses the one-form / = lady" as a Pfaff form 

I = WkdZ'' (4.2) 

of six scalar fields Z'^ and Wk ■ According to Pfaff 's theorem ||], four scalar potentials A,B,C,D are sufficient to 
describe an arbitrary covector in a four-dimensional space: 

L = AB^a + CD^a ■ (4.3) 

However, the representation of la by six potentials Wk, is more useful because it has a clear physical interpretation 

The null dust action 

5^°[Z^W^,;7a/5]= j d^L^'^iy) (4.4) 
is a functional of our six state variables, and of the metric ja/S ■ The Lagrangian density is taken in the form 

i^° = -^l7r/V%^^, (4.5) 

where la is an abbreviation for the expression (4.1). 

The equations of motion follow from the variation of the action with respect to Wk and Z'' : 

rcND 



ND 



Sz^ -{h\'^'W,l'^),c- (4.7) 
By multiplying Eq. (4.6) by Wk , we learn that ^" is a null vector field: 

l^la^Q. (4.8) 
Equations (4.6) reassert that Z'' are comoving coordinates. Equation (4.7) tells us that the three currents 

Jk = Wk r (4.9) 

satisfy the continuity equations 

VaJ^ =VaiWkn=0. (4.10) 

Because each of the three covectors Z'^ ^a is perpendicular to the null vector P, none of them can be timelike. 
If only one of the coefficients Wk in the decomposition (4.1) of P does not vanish, the congruence is hypersurface 
orthogonal (cf. Eq. (3.28)) and thus nontwisting. The covector Z'^^a is then null. In the general case of a twisting 
congruence, all three covectors Z'^^a must be spacelike: If any Z'^ ^a , k fixed, were null in an open neighborhood U, 
then I'^Z^ ^a — would imply la oc Z^ ^a in i^, and the congruence would not be twisting in U. Some covectors Z'^ q, 
can possibly become null only in lower-dimensional {d = 0, 1,2) regions of Ai. A covector Z'^^a , k fixed, can also 
become null on a 3-dimensional null hypersurface Z'^ — const on which the two remaining coefficients Wi , i k , 
vanish: Wi — . Then, of course, la oc Z'^^q, simultaneously lies in this hypersurface and is orthogonal to it. 

In Appendix C, we give two examples of twisting null congruences (one of them is the familiar ingoing principal 
null congruence in the Kerr spacetime), and illustrate the decomposition (4.1) of their tangent null covectors. The 
spacelike character of the covectors Z^^a is exhibited everywhere except in regions where the twist vanishes. 

It now becomes understandable why it would not be useful to represent la by more than six potentials. If, say, we 
wrote la — WgZ'^^a , s = 1,2,3,4, one of the spatial vectors Z" ^a could always be written as a linear combination of 
the remaining three vectors Z'' ^a , k — 1,2,3, and the decomposition (4.1) would be regained. 

The variation of the action (4.4), (4.5) with respect to the metric jap yields the energy-momentum tensor 

T"^ = 2171-1/2 (55^D/57„^. (4.11) 
Because P is null. Eq. (4.8), this tensor has the structure (2.13). The Pfaff form (4.1) satisfies the identity 
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^0{U^) = -WkAZ^pl") + Z\^^p{Wkl^) + \^c.{hl^) ■ (4.12) 

The equations of motion (4.6)"(4.8) then imply the energy-momentum conservation law. In fact, it is well known that 
the energy-momentum conservation follows from the equations of motion because of the invariance of the action (4.4), 

(4.5) under spacetime diffeomorphisms (see, e.g., jl)). We have already seen that the energy-momentum conservation 
implies that the particles of the null dust move along geodesies, Eq. (2.14). This demonstrates that our action 
(4.4)-(4.5) correctly reproduces the motion of the null dust on a given background (Al,7) . 

The null dust is coupled to gravity by adding its action to the Hilbert action 

^°[7a/3]= / d'yi°, (4.13) 
Jm 

i« = \l\"'R{y; 7] (4.14) 

constructed from the curvature scalar R{y\ 7] The variation of the total action S — + with respect to the 
metric 70/3 yields the Einstein law of gravitation (2.3) with the null dust source (4.11). The conservation law (2.4) 
then follows independently of the equations of motion directly from Eqs. (2.3) through the Bianchi identities. 

V. SPECIAL PARAMETRIZATIONS AND NULL DUST ACTION 

The geodesic equation (2.14) which follows from the action (4.4)-(4.5) is not given in affine parametrization. Rather, 
the vector field l'^{y) is chosen such that it absorbs the mass distribution M of the dust and leads thereby to the 
energy-momentum tensor (2.13). Let us now show that from any solution Wk{y) and Z^{y) of the Euler equations 

(4.6) -(4.7) of the action (4.4)-(4.5) one can construct a vector field k°'{y) given in generic affine parametrization. 
Start on a spacelike hypersurface S transverse to the dust lines P. Parametrize E by the dust space coordinates 

of points z e 5. As long as there is any dust on S , Wk{z) cannot be a zero covector in T*S. Choose an arbitrary 
vector field A^{z) £ TS such that A'^{z)Wk{z) > 0. Evolve the fields Z''{y),Wk{y) from their initial values z'^ and 
Wk{z) on S by the Euler equations (4.6)-(4.7) and define 

A{y):^{A'^{Ziy))W,{y))-\ (5.1) 

The Euler equations imply that 

V„(A-i/")=0. (5.2) 

By comparing Eq. (5.2) with Eq. (3.11), one sees that A{y) is a scaling factor (3.10) which takes into an affinely 
parametrized k" . We already know, Eq. (3.23), that the most general scaling factor A{y) can differ from our particular 
scaling factor A{y) only by a multiplicative function Aq{z) of comoving coordinates: 

A(y)-Ao(Z(y))A(y). (5.3) 

Equations (5.1), (5.3) specify an algebraic procedure by which, from any solution Z^{y) , Wk{y) of the Euler 
equations (4.6)-(4.7), one can construct the most general scaling factor A{y) which takes the covector field la ~ 
WkZ'^^ct into a covector field 

fc„(2/)=A(y)/„(y) (5.4) 

in affine parametrization. Equation (5.4) simultaneously tells us how to scale the potentials Wk into the corresponding 
potentials Wk of the Pfaff form of fc^ : 

ka — WkZ^ _a , with Wk — AWk ■ (5-5) 



■^The mixed brackets in R{y; 7] indicate that the curvature scalar 7? is a function of y and a functional of ja/jiy') ■ This 
convention is used throughout the paper. 
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The mass distribution 



M = (5.6) 

associated with the affine parametrization (5.4) satisfies a continuity equation (3.17). The potentials Wk associated 
with la also satisfy the continuity equation (4.10). However, because in general V^^" 7^ 0, the potentials Wk do not 
stay the same along the dust lines: 

I'^VaWk^O. (5.7) 

On the other hand, by virtue of the continuity equations (4.10) and (3.17), the potentials Wk associated with an 
affinely parametrized of Eq. (5.4) do stay the same along the dust lines: 

A:"VaWfc = 0. (5.8) 

Equations (3.7), (3.9) and (5.5) enable us to interpret Wk geometrically as projections of the null field ka into the 
hypersurfaces — T^pp(w, z), v — const, of afRne foliation: 

v.k^ka'-^3f^, (5.9) 

Notice that while ka in afhne parametrization is built from a solution (jj) , Wk (y) of the Euler equations by 
differentiations (5.5) and algebraic manipulations (5.1), (5.3), the construction of the affine parameter v — V{y) itself 
requires solving a differential equation k°'Va — 1 , i.e., an integration (3.19). 

The other special parameters, the luminosity distance (3.33) and the parallax distance (3.34), can be obtained from 
Z'^{y) , Wk{y) by algebraic operations and differentiations. The luminosity distance A(y) is simply the scaling factor 
(5.3). The parallax distance p is the reciprocal value of 



e 



A(V„P) = A\j\-'/^ (|7r/' 7"^WkZ\p) . (5.10) 



So far, we have shown how to construct the covector field ka in affine parametrization from a solution Z^(y) , Wk{y) 
of the Euler equations (4.6)-(4.7) of the action principle (4.4)-(4.5) written in the la parametrization. Let us now 
show how to enforce affine parametrization directly from an action principle. Require one of the potentials Wk in 
the action (4.4)-(4.5), say M-.— Ws, to be positive, and drop the index from the associated comoving coordinate: 
Z : = . Introduce wa ■— Wa/W^ in place of the remaining two potentials Wa , A ~ 1,2, and write the Lagrangian 
(4.5) in terms of the new variables M, Z, Z^, wa '■ 

L^° = -i|7r/'M7"'^fcafc/3, (5.11) 

with 

fc„:=Z„+WA^^,a. (5.12) 

The Pfaff form corresponding to ka is now constructed only from five potentials Z, Z^ , wa , though the action (4.4), 
(5.11)-(5.12) still depends on six scalar variables, due to the presence of M in the Lagrangian (5.11). By varying the 
action with respect to Z one obtains the continuity equation (3.17). By using the other field equations, one easily 
derives Eq. (3.16) for affinely parametrized geodesies. 

The Lagrangian (4.5) is special because it leads to the simplified form of the energy-momentum tensor, while the 
Lagrangian (5.11) is special because it leads to an affinely parametrized k". By building an additional redundancy into 
the Lagrangian, one can reach the generic form (2.1), (2.9) of the energy-momentum source. One simply introduces 
the seventh scalar M while keeping [/„ as the Pfaff form of six scalar fields Z'', Wfe : ^ 

L""^ = ~^h\'/^Mr^UaUp, (5.13) 



^By comparing Eqs. (5.13)-(5.14) with Eqs. (4.2) and (4.5), we see that Wk = M^^^Wk ■ 



9 



with 

C/„:=WfeZ^,. (5.14) 

The new Lagrangian density and all equations of motion arc then invariant under the gauge transformation (2.10), 
(2.11), 

Wfe^Wfe=AWfe (5.15) 

and 

M->M = A"2m, (5.16) 

where A{y) > is an arbitrary scaling factor. 

The canonical form of the action is the same whether one starts from the original Lagrangian (4.5), (4.1) or the 
redundant Lagrangian (5.13)-(5.14). The canonical variables recombine the redundant potentials in such a way that 
the information about the split of la into M and Ua gets lost: From the canonical variables one can reconstruct 
only la ■ It is thus not worth the effort to complicate the spacetime Lagrangian by striving to achieve a superfluous 
generality. Having learned this lesson, we take the spacetime action (4.4)"(4.5) with la given by Eq. (4.1) as our 
starting point. 

VI. CANONICAL DESCRIPTION OF NULL DUST 

The familiar ADM algorithm for casting a covariant action into Hamiltonian form works for the null dust in a 
similar way as for the ordinary dust of massive particles [Q . One foliates the spacetime M by spacelike hypersurfaces 
S, 

Y:RxT,->M by {t,x) ^ y ^Y{t,x) . (6.1) 

In local coordinates x^, a = 1, 2, 3 on E and y", a = 0, 1, 2, 3 on M, the fohation is represented by 

(t,a;")^y"==r"(t,a;'^). (6.2) 

A transition from one leaf E of the foliation to another is described by the deformation vector Y" := dY°' / dt . Its 
decomposition into the normal n" and tangential Y" directions to the leaves yields the lapse function N-^ and the 
shift vector N"-; 

Y" = N^n"^ + N^Y^^a- (6.3) 
On each leaf, the spacetime metric ^apiv) induces the intrinsic metric 

gab{t, x) = -iap{Y{t, x)) y",,(i, x)YP^i{t, x) . (6.4) 
The spacetime metric is reconstructed as 

7"^ = -n"n'3 + , (6.5) 

where g""^ is the inverse of gab , and the determinants I7I of jap and |(?| of gab are related by 

|^|l/2^^^|^|l/2^ (6.6) 

Scalar fields on Ai, such as the null-dust variables Z'', Wk can be pulled back to JR x E by the mapping (6.1). By 
using Eq. (6.3) we obtain 

= (7V^)^V\ (6.7) 

where we have introduced the normal velocities 

:= - .aN" , Z^ ^a ^ ^aY" .a ■ (6.8) 
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This allows us to write the null-dust action (4.4)-(4.5), with 1^ given by (4.1), as an integral over 1? x S, i.e., in the 

(3+l)-split form: 

5N°[Z^W^fc; 5„6,iV^,7V«] = [ dt [ cfxl""^. (6.9) 

JM JS 

The Lagrangian density L^^ on JR x S is a quadratic form of the Lagrange multipliers Wk ■ 

iND = 1 |^|l/2 (^(^^±ylyiyj _ ^±gij^ _ (g_;^0) 

The metric 

g'J{t,x)=g'^^Z\aZ\b (6.11) 

is the induced metric on S expressed in the basis of comoving coordinates Z*. 

By varying the action with respect to Wj , we get a system of linear homogeneous equations for Wj : 

{g'^ - {N^y'^V'V^) Wj=Q. (6.12) 

This has a nontrivial solution only if the determinant 

det {g'^ - {N^y'^VW^) = (l - {N^y'^gi^VW^) det [g'^) (6.13) 

vanishes. This imposes the constraint 

gi^V^V^ - {N^f = (6.14) 

on the velocities Z'- . (Here, gij is the inverse of g^^ . We can use it for lowering the dust space indices.) 

If the constraint (6.14) is satisfied, Eq. (6.12) has a solution Wj oc V,- . Of course, the homogeneous equation (6.12) 
determines only the direction of Wj , leaving W = g^^WiWj undetermined. We write the general solution in the form 

Wj = VWVj {ViV')-^/^ , (6.15) 

where W is an arbitrary positive factor. 

By substituting this solution (6.15) back into the Lagrangian (6.10), we eliminate from the action the multipliers 
Wj , replacing them by a single multiplier W : 

^ND ^ 1 |^|i/2 ^ ({N^y^gijVW^ - N^) . (6.16) 

The reduced action is a functional of W and Z'^. Its variation with respect to W reproduces the constraint (6.14) 
which enabled us to express Wj in terms of W and Z^ , Z^ , Eq. (6.15). Its variation with respect to Z'^ gives an 
equation which, modulo the constraint (6.14) and Eq. (6.15) considered as a definition of Wj , is equivalent to the 
equations of motion obtained by varying the original action (6.9)-(6.10) with respect to Z'^. The reduced action 

S^^[Z'',W;gab,N^,N'']=[ dt ( d^xL^^ (6.17) 

J m, Jt, 

with the Lagrangian (6.16) is thus entirely equivalent to the original action (6.9) with the Lagrangian (6.10). 

In order to bring the reduced action to canonical form, wc perform the Legendre dual transformation from (Z*^, Z*^) 
to (Z'^,Pfc), leaving W as, a, multiplier. First, we introduce the momenta 

QrND 

Pk ■■= = \g\'/'W{N^)-'Vk . (6.18) 

To clarify their physical meaning, we return to the definition (6.15) of Wj , the decomposition (4.1) of , and Eqs. 

(6.7)-(6.8) for the normal velocity. In this way we learn that Pk are normal projections of the currents introduced 
inEq. (4.9): 

Pk = \9\^^^J^na. (6.19) 
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Equation (6.18) can be inverted to obtain the velocities 



= N^\g\-^''^W-^g^^Pj + N'^Z^^a ■ (6.20) 

This leads to the Hamiltonian 

i/ND _ p^^fc _ ^ND ^ N^Hf + N'^H^'' (6.21) 
which is a linear combination of the momentum density 

H^"" = PkZ\a (6.22) 



and the energy density 



hT = \w-^\gr^'^g''P.Pj + \w\g\"^ (6.23) 

= \w-'W'g'''HrHr + \w\g\"' (6.24) 
of the dust. The canonical form of the action then reads 

S''''[Z'',Pk,W;gab,N^,N'']^ J^dt J^d'x {^P^Z'^ - N^Hf"" - N'^H^'') , (6.25) 

where H^^ and Hf^ are given by Eqs. (6.22) and (6.24). 

At this stage, we are finally able to eliminate the last remaining multiplier W. By varying the action (6.22)-(6.25) 
with respect to W, we obtain an equation 



SW dW 
which determines W in terms of the canonical data: 



^"^=0 (6.26) 



W = \g\-"^ - \g\-'^^ ^g'^'H^H^'' ■ (6.27) 

By substituting this solution back into we obtain 



Hr = ^/g'^'H^^H^'' ■ (6.28) 

We see that W is just the scalar form W — \g\^^^^H^^ of the Hamiltonian density . The final expressions 
(6.22) for the momentum density and (6.28) for the energy density are simple: The form of the momentum density is 
dictated by the requirement that it generate the Lie derivative change of the scalars Z'^ (x) and scalar densities Pk (x) 
under spatial diffeomorphisms LDiffS The energy density is the norm of the momentum density with respect 

to the spatial metric. The resulting reduced canonical action S^^[Z'', Pk] gab, N-^^N"-], with (6.22), (6.28) yields the 
Hamilton equations for Z''{t,x) and Pk{t,x). These describe the evolution of the null dust on a given geometrical 
background 7a/3 ^ {N-^{t,x), N°'{t,x), gab{t,x)). ^ 

From the solution of the Hamilton equations we can reconstruct the null vector I" which provides the spacetime 
description of the dust. It holds that 

r = l^n°' + rV" , (6.29) 
where l'^ and la are expressed as functions of the canonical variables: 



*Our Lagrangian and Hamiltonian formalism can easily be generalized to several mutually noninteracting species (streams) 
of null dust. This may be useful for the canonical treatment of spherical collapse, in which the ingoing null dust is turned into 
an outgoing null dust at the center of symmetry, or for the canonical treatment of models involving colliding streams of dust 
with plane or cylindrical symmetry (see references in [Q). 
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/-L = -l^n" = (6.30) 

/, = ;„y",, = M/-i/2|5r^/X''°- (6.31) 

Here, of course, W stands for the scalar form (6.27) of the energy-momentum density. One can check that I" is a 
null vector by virtue of its construction (6.29)-(6.31): 

l^'l^^ -{l^)^ + n^ = 0. (6.32) 

The background variables N^{t,x), N°'{t,x) and gab{t,x) in the dust action S^^[Z'' , Pk ; gab ,N^,N°-] are not to 
be varied. The Hamiltonian formalism for null dust on a given background is thus entirely unconstrained. To couple 
null dust to geometry, we must add its action S^^ to the gravitational Dirac-ADM action 

5^ [5ah,p"'; N^,N-] = [ dt [ d?x {p-'gab - N^Hf - N^H^) (6.33) 

with the standard gravitational super-Hamiltonian and supermomentum densities 

Hl{x- gab,p'''] = Gabcd{x-g) p'^\x)p-\x) - \g\"^R{x- g] , (6.34) 

Gabcd = 2 1^1"^^^'^^'''=^'"^ Sadffbc ~ 9ab9cd) , (6.35) 

H^ix; gab^p""'] - -2Dbpl{x) ; (6.36) 

here, Db is the spatial covariant derivative. 

The variation of the total action with respect to the lapse and the shift N°- then leads to the familiar Hamiltonian 
and momentum constraints 

Hi_ := Hf + = , (6.37) 
Ha + H^^ = (6.38) 

for the coupled system. 

VII. NULL DUST CONSTRAINTS THAT GENERATE A LIE ALGEBRA 

By using the supermomentum constraint, one can replace the momentum density of the dust by the gravita- 
tional density in the expression (6.28) for the dust energy density . This brings the constraint system (6.38), 
(6.37) into an equivalent form (6.38) and 



Hi_ := Hf + ^g-<'H2H^ - . (7.1) 

Only the supermomentum constraint (6.38) contains the dust variables. The new Hamiltonian constraint (7.1) is 
constructed solely from the gravitational variables gab , p""^ ■ Alternatively, one can get rid of an inconvenient square 
root by rewriting Eq. (7.1) in the form 

G := {Hff - g'^'^H^H^ - . (7.2) 

Under the positivity condition 

- H^^ = H^<0, (7.3) 

the constraint (7.2) is equivalent to the constraint (7.1). 

Brown and Kuchai' ||l| proved a remarkable fact that the densities (7.2) have strongly vanishing Poisson brackets: 

{G(x), G(x')} = 0. (7.4) 

By coupling gravity to other simple sources, Kuchaf and Romano [O and Brown and Marolf |l2) produced other 
densitized expressions constructed from the scalar variables g~^^^ HY and g~^g°''^H^H^ which also have strongly 
vanishing Poisson brackets. Markopoulou psl posed the question what is the most general density 



13 



of weight w constructed from these variables which has the strongly vanishing Poisson brackets 

{Fix), F{x')}^0. (7.6) 

She found an algorithm for generating all such densities. The density (7.2) of weight 2 still seems to be the simplest. 
Among others, there is the scalar form 



G, Hf + Jg^^H2H^ = (7.7) 



of the constraint (7.1) which, as we have just seen, describes null dust. 

The constraints H2 = = of vacuum gravity can be replaced by an alternative system 

ijG = = G (or i/„G = = G^) . (7.8) 

Unlike the original constraints, and G (or and G ^ ) generate a true Lie algebra. Unfortunately, in vacuum 

gravity the new constraints (7.2) (and similarly (7.7)) do not generate the evolution of the geometric data gab , p""^ 
into a Ricci-flat spacetime. Expression (7.2) is flawed because its Hamiltonian vector field vanishes on the constraint 
surface (7.8), while the expression (7.7) is flawed because its Hamiltonian vector field is ill-defined for — . 

No such difficulty exists for null dust. The momentum constraint (6.38) is different from the vacuum constraint 
= and, as long as there is any dust at the point in question, H^^ and hence cannot vanish. The Hamiltonian 
vector fields of the dynamical variables (7.2) or (7.7) then do not vanish on the constraint surface (6.37)-(6.38) of 
the null dust coupled to geometry. The new constraints (7.2) or (7.7) correctly generate the evolution of geometry 
produced by null dust. Moreover, as in vacuum spacetime, the constraints (6.38) and (7.2), or (6.38) and (7.7), 
generate a true Lie algebra. It is thus advantageous to bring the constraints to one of these forms before attempting 
to quantize the coupled system. 

Why is it that the presence of the null dust does not affect Eqs. (7.2) and (7.7) that hold in vacuum gravity? The 
energy-momentum tensor of the null dust satisfies the condition 

j,a^j,7/5^0. (7.9) 

Conversely, any symmetric tensor T"^ which satisfies Eq. (7.9) must either vanish, or there exists a null vector P 
such that 

T"^ = (7.10) 

The Einstein law of gravitation (2.3) then implies that P is a geodesic vector field, i.e., the Euler equations of motion 
for the null dust. The simple tensor equation 

QO^^Ql/S ^ (7.11) 

imposed on the Einstein tensor thus ensures that the geometry 7q,^ is necessarily produced by null dust according to 
Einstein's law of gravitation.^ 

The _L_L projection of Eq. (7.11) gives 

{G^^f-g^'G^aG^b^Q. (7.12) 

Because the _L_L and _L|1 projections of the Einstein tensor yield the gravitational super-Hamiltonian and supermo- 
mentum [flOi, 



^Equation (7.11) is perhaps the simplest example of the Rainich-type geometrization of a source field. The general task is to 
find equations for the Einstein tensor which are equivalent to the Einstein law of gravitation together with the field equations 
for a given source. The problem was first formulated for the Einstein-Maxwell system by Rainich and solved by him under 
the assumption that the electromagnetic field is not algebraically special (null) jl4|, ||l^. The Rainich problem for the null 
electromagnetic field was solved by Hlavaty ||l^. The much simpler scalar field case was analyzed by Peres and by Kuchaf 
[|l8|. The spinor field was treated by Kuchai^ |19| and the Proca field by Bicak pojl. 
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= -\9-"'lil , = \9-'l^ll^^ , (7.13) 

Eq. (7.12) is equivalent to the constraint (7.2). We have already noticed that under the energy positivity condition 
(7.3) the constraint (7.2) is equivalent to the constraint (7.7). The Rainich-type condition (7.11) thus connects the 
new form (7.2) or (7.7) of the Hamiltonian constraint with the spacetime picture. 



VIII. CONSTRAINT QUANTIZATION OF GEOMETRY COUPLED TO NULL DUST 

We have cast the constraint system for geometry coupled to null dust into a form in which it generates a Lie algebra. 
In this process, the Hamiltonian constraint has been replaced either by the constraint (7.2) or by the constraint (7.7). 
Either of these constraints have vanishing Poisson brackets (7.4) . The momentum constraint is left in its original form 
(6.38) and (6.22): 

-ffa(x) Pfe(.T)Z\,(.T) + ijG(a;) = . (8.1) 

The momentum constraints (8.1) close in the way characteristic for the Lie algebra LDiffS of the diffeoniorphisni 
group DiffS. The Poisson brackets of G[x) (or G ^ix)) with Ha{x') close into G{x) (or G ^{x)) in the way which 
reflects the transformation behavior of G(x) ( or G ^ (x) ) under spatial diffeomorphisms DiffS : G{x) is a density of 
weight 2, while Gy (x) is a scalar. 

As for ordinary dust, the constraint system can be vastly simplified by the introduction of an alternative set of 
canonical variables which reflect the fact that the dust particles define a preferred system of coordinates on S. The 
mapping Z : S — > 5 which, in local coordinates, assumes the form 

zk^Z'^ix"), (8.2) 

takes the tensorial variables gab{x) and p°'^{x) on S into corresponding tensors gij{z) and p^^ (z) on the dust space S : 

g,^{z) ■.= X\,{z)X\,{z)gab{X{z)) , (8.3) 
dX{z) 



Z\aiXiz))Z\,{Xiz))p''''iXiz)) . (8.4) 



dz 

Here, the f-dependent mapping X : S ^ is simply the inverse of Z, 

X := Z-\ (8.5) 

and \dX(z)/dz\ is the Jacobian for the change of variables 2;" — X'^(z). 
We rewrite the supermomentum constraint (8.1) in the form 

H^kix) ■.^Ha{x)Zl{x) = Pk{x)+H^{x)Zl{x) = 0. (8.6) 

Here, 

Zl{x):=X\k{Z{x)) (8.7) 

is the inverse matrix to Z^ ^a{x) : 

Z\,{x)Zl{x)^5i. (8.8) 
The new supermomentum H-\k{x) smeared by a new shift N^^{x), 

i7T[iVT]:= [ £{x)N^\x)H^k{x), (8.9) 

generates through the Poisson bracket the change 

Z'^ix) := {Z^ix), H^[N^'j = N^'^ix) (8.10) 

of the dust coordinates Z''{x) by the amount N^'^{x) 
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One can prove that the 5-variables g^{z), p^^{z) along with the dust frame variables Z^{x) and the new supermo- 
mentum H-^k{x) form a canonical chart [y . In particular, this means that the new constraint functions Pfc(a;) := H-^k{x) 
have vanishing Poisson brackets among themselves and are the momenta Pfe(a;) canonically conjugate to the dust frame 
variables Z^{x). Further, because the Poisson brackets of the 5-tensors g^j{z), p^^ (z) with the smeared supermomen- 
tum (8.9) vanish, these 5-tensors are invariant under the shifts (8.6). 

In terms of the new canonical variables g^j{z), p^^ (z) and Z^{x), Pfc(x) the momentum constraint (8.6) reduces to 
the condition that the canonical momentum Pfc(x) vanishes: 

Pk{x) = 0. (8.11) 

The Hamiltonian constraints (7.2) or (7.7) can then be mapped to the dust space S according to their weight: 



G(z): 



dX{z) ' 



G(X(z)) = 0, (8.12) 



dz 

G^(z):=G^(X(z)) = 0. (8.13) 

The 5-constraints (8.12)-(8.13) are the same functionals of the 5-tensors gij{z) and p^-' {z) as the E-constraints 
(7.2) or (7.7) were of the S-tensors gab{x) and p°'''{x) . In other words, G(z) is obtained from G{x) and G y(z) is 
obtained from G ^{x) by replacing the S-tensors gab[x), p"-^{x) by the corresponding 5-tensors gij{z), p^^{z). The 
5-constraints (8.12) and (8.13) have strongly vanishing Poisson brackets, i.e., they generate an Abelian algebra. 

In the Dirac method of quantization, constraints are turned into operators and imposed as restrictions on the state 
functionals of the system. We choose to work with the dust space variables, so the quantum states of the system are 
functionals ^[Z^ g] of the canonical coordinates Z^{x) and g^^jiz), and the constraint operators Pk{x) and G{z) (or 
Gy(z)) are obtaned by quantizing the classical expressions (8.11)-(8.13). The transition is easy for the momenta 
(8.11) which are simply replaced by the variational derivatives 

P^(-) - • (8.14) 

The operators (8.14) automatically commute, 

P,(a;), Pj(x')] =0. (8.15) 

It is far from clear how to replace the remaining classical constraints by operators which not only commute with 
Pfe(x), but also among themselves. We shall proceed under the assumption that there exists a factor ordering and 
regularization of G:—G{z; g^j{z), p^-' {z)] and/or G^:=G^{z; g^j{z), p'''' {z)] which achieves this goal. If so, the 
constraint operators can consistently annihilate the physical states. The momentum constraint 

Pk{x)^[Z,g]^0, (8.16) 

where Pfc(a;) is interpreted as the variational derivative (8.14), means that the state functional g] cannot depend 
on Zf^ix): 

* = • (8.17) 

The constraint system is thereby reduced to a single oo"^ nontrivial condition that G := G(z; (z), p'-' (z)], 
(or G y :=G y(z; g^j{z), p*"'(z)] ) annihilates the state functional: 

G(z;g,^-(z),r''(z)]*[g] =0. (8.18) 

The caveats which need to be born in mind when implementing such a formal procedure for gravity coupled to ordinary 
dust are carefully spelled out in [Q . An additional difficulty with null dust is that there is no natural variable which 
would play the role of internal time. As a result, unlike for ordinary dust, the quantum constraint (8.18) does not 
have the form of a functional Schrodinger equation. It is thus unclear how, even formally, to turn the space of its 
solutions into a Hilbert space. 



16 



IX. COMPARING NULL DUST WITH ORDINARY DUST 



Ordinary dust coupled to gravity was turned into a Hamiltonian system and formally quantized by Brown and 
Kuchaf This scheme turns out to be both similar to and characteristically different from the description of null 
dust given in this paper. We shall outline the basic similarities and emphasize the differences. 

The spacetime action 

5° [T, ; M, Wk ; 7a/3] ^ I d^yL^{y) (9.1) 

J M 

of ordinary dust is constructed from eight scalar fields , Wk and T, M. The Lagrangian density L^{y) has the 
form 

i°--^l7l'/'M(7"%„C/^ + l). (9.2) 
The four-velocity is expressed as the Pfaff form 

C/a - -r„ + WkZ\o. (9.3) 

of seven scalar fields W/c, and T. The matter equations of motion are obtained by varying the dust action 
(9.1)-(9.3) with respect to the state variables M, Wk: T and Z'': 







5S^ 

5M ~ 2 



l\l\'^'{r''U^Up + l), (9.4) 



= ^ =-|7|i/2MZ^„{/", (9.5) 



SWk 
6T 
SZ^ 



(|7|i/2Af{/"),„, (9.6) 
|7|i/2MVFfei7"),„. (9.7) 



They lead to the interpretation of the state variables. Equation (9.5) is analogous to Eq. (4.6) for null dust. It ensures 
that the three vector fields Z*^ are constant along the fiow lines of C/" and therefore their values z*^ can be interpreted 
as comoving coordinates for the dust. Equation (9.4) ensures that the four-velocity t/" is a unit timelike vector field. 
It is analogous to Eq. (4.8) which guarantees that the four-velocity P of null dust is lightlike. Equation (9.6) allows us 
to interpret M as the rest mass density of the dust and expresses the law of mass conservation. It is analogous to Eq. 
(3.17) for the null dust in affine parametrization. Equation (9.7) can be interpreted as the momentum conservation 
law. It is analogous to Eq. (4.10) for the null dust written again in affine parametrization. By multiplying Eq. (9.3) 
by U" and using the field equations (9.4)-(9.5), we learn that 

r„c/" = 1 , (9.8) 

i.e., that T is the proper time between a fiducial hypersurface T = and an arbitrary hypersurface T = const along 
the dust worldlines. From Eq. (9.3) we see that the Wk variables are the projections of the four-velocity Ua to the 
hypersurfaces of constant T expressed in the dust space cobasis Z^ ■ Due to the conservation laws (9.6)-(9.7), these 
projections remain the same along a fiow line of U" . In comparison, Wk for the null dust is the component of the 
null covector P in the dust space cobasis Z^ ■ These components are not conserved along the fiow lines, Eq. (5.7). 
However, when one rescales I" into an affinely parametrized k" by Eqs. (5.4)-(5.5) and projects into hypersurfaces 
of constant affine parameter one obtains the components Wk of Eq. (5.9) which are conserved along the fiow lines, 
Eq. (5.8). 

The main difference between the actions and is that the dust action depends on eight variables T, M and 
Z^ : Wk , while the null dust action depends only on six variables Z^ , Wk ■ The interpretation of the variables Z^ as 
the comoving coordinates and Wk as the projections of the four-velocities (or P) into hypersurfaces of constant 
T (or U) is analogous. The variables T and M do not appear in the null dust action (4.1), (4.4)-(4.5). This refiects 
the fact that the mass function M of the null dust is not uniquely determined and it was absorbed into the definition 
of Z". Similarly, the affine parameter along the null geodesies is not uniquely determined. If one chooses to enforce 
the affine parametrization by taking the null dust Lagrangian in the form (5.11), the corresponding M occurs in the 
action, but the Pfaff form of an affinely parametrized k^ , Eq. (5.12), contains only two independent scalars w^ . One 
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can work in a totally arbitrary paramctrization by letting the Lagrangian density to depend on seven variables M, , 
Wft instead of six, Eqs. (5.13)-(5.14), but then the action becomes gauge invariant under the scalings (5.15)-(5.16), 
which makes it effectively dependent only on six of these variables. 

These similarities and differences are reflected in the canonical form of the action. For ordinary dust, the energy 
density and momentum density depend on four pairs of canonical variables, T and P, and Z'^, P/. . They 
take the form 

= PTa + PkZ^^^ (9.9) 

and 



= \/P^ + g'^bR^Hf . (9.10) 

On the other hand, similar expressions for null dust, Eqs. (6.22) and (6.28), depend only on three pairs of canonical 
variables, Z^ and Pk ■ This difference is vital. While ordinary dust has four degrees of freedom per space point a; G E, 
null dust has only three. 

The rest mass density M of ordinary dust is directly related to the momentum P : 



The mass function and affine paramctrization of null dust arc ambiguous and their only invariant c;ombination is the 
null vector I". This can be reconstructed from the canonical data, Eqs. (6.29)-(6.31), rather than the mass function 
and the four-velocity separately. 

Formally, the momentum and energy densities (6.22) and (6.28) of the null dust arc obtained from the corresponding 
expressions (9.9)-(9.10) for ordinary dust simply by putting P = and forgetting all about its conjugate variable T. 
This should not hide the fundamentally different ways in which the Dirac-ADM action is obtained from the spacetime 
action. The Lagrangian (9.2)-(9.3) for ordinary dust is nondcgcnerate in the velocities T, Z'^. The expressions for 
the momenta P, Pk can be inverted to yield the velocities. The momenta are in a one-to-one correspondence with 
the multipliers M and Wk and hence their variation yields equivalent equations. The spacetime action, so to speak, 
is in an 'already parametrized form'. 

To cast the spacetime action (4.4)-(4.5) of the null dust into canonical form requires an entirely different procedure. 
The null dust Lagrangian (6.8), (6.10) is singular in the velocities Z''. The definition equations for the momenta Pk 
cannot be inverted. They yield three constraints 

6'^''WjPk = (9.12) 

demanding that the multipliers Wk be parallel to the momenta Pk , which leaves the magnitude of Wk undetermined. 
The variation of the action with respect to Wk leads to the constraint (6.14) on the velocities Z''. If this constraint is 
satisfied, the multipliers Wk can be replaced by a single multiplier W and the Lagrangian L^^ cast into an equivalent 
form (6.16) which is regular in the velocities. This allows one to perform the Legendre dual transformation to the 
canonical form of the action. The final elimination of the multiplier W (analogous to the final elimination of the 
mass multiplier M from the canonical action for ordinary dust) leads to the null dust momentum and energy densities 
(6.22) and (6.28). To summarize, though these final expressions have similar structure as the densities (9.9)-(9.10) 
for ordinary dust from which they can be obtained by putting P = 0, their derivation is fundamentally different. 

After the dust is coupled to geometry, the parallels and differences between ordinary and null dust are brought into 
a new perspective. The momentum and Hamiltonian constraints for ordinary dust can be resolved with respect to 
the four dust momenta P, Pk which brings them to an equivalent form 

H^k:^Pk + Z^H^ + VGTaZ^ = 0, (9.13) 
H^:=P-VG^0, (9.14) 

where G is given by Eq. (7.2). The new constraint functions H^k = {H^ ,H'\k) have strongly vanishing Poisson 
brackets: 

{H^k{x),H^l{x')} = Q. (9.15) 

The imposition of the constraints (9.13) and (9.14) as operator restrictions on the states '^\T,Z^; gab^P""^] leads 
to a functional Schrodinger equation with formally conserved inner product. By mapping the constraints into the 
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dust space, the momentum constraint is eliminated and what remains is a single functional differential Schrodinger 
equation 



(P(z) - ^G{z- g, p]) ^[T{z),g{z)] = 0. (9.16) 

The nuU dust constraints in the form (8.6), (7.2) can again be obtained from the ordinary dust constraints (9.13)- 
(9.14) by disregarding the canonical pair T, P (and squaring Eq. (9.14)). By mapping them into dust space, the 
momentum constraint is again eliminated. By imposing the only remaining constraint as an operator restriction on 
quantum states, one again gets a single functional differential equation (8.18). However, and this is an important 
difference, Eq. (8.18) is not a Schrodinger equation like Eq. (9.16) because there is no internal time T(z). It is thus 
not clear how to introduce an inner product in the space of its solutions. 

Both ordinary dust and null dust provide a standard of space in canonical gravity because the dust particles 
introduce into spacetime a privileged dust frame S labeled by comoving coordinates Z''{x). The crucial difference 
is that ordinary dust provides also a standard of time: It has an additional degree of freedom T{x) which can be 
physically interpreted as the proper time along the dust worldlines. Null dust does not have any corresponding degree 
of freedom because afffne parametrization of null geodesies is ambiguous. It thus fails to provide a standard of time 
to the spacetime in which it moves. The story of ordinary dust is that of time regained. The story of null dust is that 
of time lost again. 
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APPENDIX A: NULL DUST AND GEOMETRICAL OPTICS 

If at each spacetime point all the energy is transported in one direction with the speed of light, it is appropriate to 
describe the matter by the energy-momentum tensor of null dust, 

T"'3 = Mk^kf^ . (Al) 

The energy-momentum tensor (Al) may be considered as representing an incoherent superposition of waves with 
random phases and polarizations but moving in a single direction.^ It is also called the 'geometrical-optics' or 'pure 
radiation' energy-momentum tensor. 

As an example, consider the Maxwell theory. (See, in particular, §22.5, for a detailed exposition of geometrical 
optics in curved spacetime.) If the electromagnetic waves can locally be regarded as plane waves propagating through 
spacetime of negligible curvature, one can write the electromagnetic vector potential Aa in the form 

Aa, = Re (a„e*®) . (A2) 

Here, in the first approximation, the complex amplitude aa{y) is independent of the wavelength and is slowly changing 
as a function of spacetime position y, while the scalar function & {y) is a rapidly changing phase. Following the standard 
procedure |^ , one introduces the wave vector 

ka = 0,a ; (-^3) 

the (real) scalar amplitude 

A-(^„yl")i/2 = (a"a„)i/2, (A4) 
and the (complex) unit polarization vector 



^This is different from the energy-momentum tensor of perfect fluid with the equation of state p — M/3, which represents the 
superposition of waves with random propagation directions. 
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ea= A '^aa ■ (A5) 

As a consequence of the source-free wave equation and the Lorentz gauge condition, both written in the first order 
of the geometrical optics approximation, the quantities (A3)-(A5) obey the following set of equations: 

/fc^r = , (A6) 

kl^W/ik" = , (A7) 

V„(yl2/fc") = 0, (A8) 

and 

fc"e„ = 0, k^Vpea = 0. (A9) 
From Eq. (A7) we see that the null vector k" is affinely parametrized. The electromagnetic field tensor is given by 

= 2Re (zAe^®fc[„e/3]) . (AlO) 
It represents the electromagnetic field of type N (the null field) since it satisfies the relations 

{Fa.p + iF*p)k'' = , K;3F"'5 = F^fiF*^P = , (Af f ) 

where F^^ is dual to Fap . Equations (A7) and (A8) imply the covariant conservation law for the electromagnetic 
energy-momentum tensor 

T"^ = A2/fc"fc^. (Af2) 

We see that the phenomenological null dust equations (2.9), (3.16)-(3.18) are the same as Eqs. (A6)-(A8) and 
(Af 2) of the high-frequency limit of the Maxwell theory if the null vector field ka is defined by Eq. (A3) and the 
mass distribution M is identified with the square of the scalar amplitude A : 

M^A^. (Af3) 

Null dust thus exhibits all features of the geometrical optics limit of Maxwell's theory except for the polarization 
properties. However, starting from a solution of the null dust equations one can always construct a polarization vector 
Ca such that Eqs. (A9) are also satisfied. This yields the tensor (AfO) which can be regarded as an electromagnetic 
field tensor in the geometrical optics approximation. 

The laws of geometrical optics can also be interpreted as describing photons that move along null rays with the 
flux vector which is determined by the amplitude A and the null vector A:" (see |^ for details). 

The lightlike particles need not necessarily be photons. It is quite obvious that similar conclusions can be reached 
for all zero-rest-mass fields in high-frequency limit. For example, by employing the geometrical optics form (AI) of 
the energy-momentum tensor, several authors ]2f| ] studied the gravitational collapse with escaping neutrinos. 

A somewhat special case is the gravitational field itself. Careful studies of the high-frequency limit of the gravita- 
tional radiation by Isaacson and others j2^ have shown that the energy-momentum tensor (Af 2) and the null vector 
field /c" which satisfy Eqs. (A6)-(A8) also describe the behavior of high-frequency gravitational waves. The metric 
tensor perturbations representing high-frequency waves are given by 

hap = Re (^(aa/3 - \a 7a/3)e*®) , 

7 a^p , (AI4) 



a 



where lap is the background metric (the source of which may be the high-frequency waves themselves) . By applying 
the geometrical optics approximation to the perturbed Einstein's equations, one arrives again at the equations (A3), 
(A6)-(A8), and (AI2). Instead of the scalar amplitude (A4) one now gets 

A^{\a^Paapf'\ (AI5) 

One also obtains the equations for the polarization tensor eap — aap/A, analogous to Eqs. (A9) (see (22|). The 
Riemann tensor of the metric (AI4) has the Petrov type N. The gravitational field in the high-frequency limit is 
null, similarly as the electromagnetic field. The well-known peeling-off property of exact radiative (zero-rest-mass) 
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fields in asymptotically flat spacetimes ^ implies that at large distances from the source these fields are null, having 
the structure of plane waves. In asymptotic regions one can even describe exact solutions of the field equations in 
terms of null dust. In such situations, one can usually find a natural parametrization of null rays - for example, by 
the proper time of distant observers at rest with respect to an isolated source. 

The variational approach of MacCallum and Taub |2^] to the high-frequency gravitational waves is especially rel- 
evant for the present paper. By applying the 'averaged Lagrangian technique' of Witham to the second variation 
Lagrangian for the perturbations of vacuum gravitational field, these authors give a variational principle for approx- 
imately periodic gravitational wave described by metric perturbation of the form (A14). Their principle, derived by 
perturbing and averaging the Hilbert action, implies the geometrical optics equations (A3), (A6)-(A8), and (A12), 
with A given by Eq. (A15). This principle is closely related to our variational principle for null dust (given in Eqs. 
(4.4), (4.5)), in the special case of the hypersurface orthogonal vector field P. 



APPENDIX B: EXACT SOLUTIONS WITH NULL DUST: EXAMPLES AND SOME RECENT 

APPLICATIONS 

As an illustration, we shall give a few examples of known exact spacetimes with null dust. (A detailed survey of 
such solutions found before 1980 is given in The cosmological solutions with null dust were recently reviewed in 
p5| , and the solution representing colliding plane gravitational waves accompanied by null dust in pGf.) 

Among the simplest solutions directly related to the fields arising in the geometrical optics limit are conformally 
flat null dust solutions representing special plane waves. They are described by the line clement (see, e.g., p^) 

ds2 = _i$2^^t_)(a,2 y^)du^_ - 2du+du^ + dx^ + dy'^ , (Bl) 

where $ is an arbitrary function of a retarded time u_ . The corresponding energy-momentum tensor is 

T^p = ^^k^kp ; (B2) 

the only nonvanishing component of the null covector ka is ku_ — 1. These solutions can always be interpreted as 
exact solutions of the Einstein-Maxwell equations with the null electromagnetic field given by Fap — 2<l>(M_)A;[Q,e^], 
where Cq. = (0, 0, cos'0, sinip) contains an arbitrary function tp = ipiu-) (cf. Eq. (AlO)). Cylindrical gravitational 
waves accompanied by null dust are also known p7[ |. 

A more complicated class of radiative solutions with 'spherical' gravitational waves and null dust is formed by 
the Robinson-Trautman solutions The energy-momentum tensor has again the form (B2), but the function <i> 
is now given by $^ = ti^(C, u_)/Tr , where C is a complex spatial coordinate, v is an affinc parameter along the 
rays, and u_ is a retarded time. The function n may be arbitrary. If, however, these solutions should represent 
exact Einstein-Maxwell fields, n must have the form n? — 2hhP^, where /i(C,(^, u_) and P{C,C,,u-) satisfy certain 
additional conditions p^ . The Robinson-Trautman solutions with null dust include Vaidya's spherically symmetric 
metric as a special case. In fact, if the evolving null dust is homogeneous, all such Robinson-Trautman spacetimes 
approach the Vaidya's metric as the retarded time goes to infinity |2^ . 

The null vector field fc" in the solutions we have mentioned is hypersurface orthogonal and the corresponding null 
congruence is thus nontwisting. The twisting null dust solutions are discussed in Q , the best known simple example 
being the 'radiating Kerr metric'. 

Some exact solutions with null dust can also be interpreted as exact solutions of Einstein's equations coupled to a 
massless scalar field ||3l|]. However, given a conserved energy-momentum tensor in the form (Al), it is not necessarily 
true that the mass distribution M and the null vector field fc" represent an electromagnetic or a massless scalar field. 
However, if the null vector field k" is shear-free, a corresponding nontrivial solution of Maxwell's equations can be 
found by virtue of the Mariot-Robinson theorem |32j . 

Recently, certain exact solutions with null dust which can be interpreted as 'relativistic rockets' have been explored 
in connection with the properties of gravitational radiation js^. A number of studies have also been devoted to 
colliding plane and cylindrical systems with null dust [ p4[ . 

Above all, as we have already stated in the Introduction, the null dust models have been recently used to clarify 
the formation of naked singularities during a spherical gravitational collapse js^], in the studies of the mass inflation 
inside black holes |^^, and in the models attempting to describe the formation and Hawking evaporation of black 
holes 
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APPENDIX C: DESCRIPTION OF TWISTING NULL CONGRUENCES BY PFAFF FORMS: TWO 

EXAMPLES 



Since null congruences are somewhat unusual, we give here two examples of twisting null congruences described by 
the scalar potentials and Wi. 

I. In a flat spacetime with Lorentzian coordinates {t, x, y, z), consider a system of lightlike particles which, in each 

plane perpendicular to the 2;— axis, move in mutually parallel straight lines. As one passes from one plane z = const 
to another, the angle a between particle trajectories and the x-axis smoothly changes with z : a = a{z) G [0, 2tt). It 
is easy to see that the null worldlines form a twisting null congruence: 

t = v + to, 

X = V cosa{z) + xq , (CI) 
y = v sma{z)+yo, 
z = Zo, 

where e 1? is an afline parameter. The tangent null vectors k" = dx^/dv are given by 

fc" = (1, cosa(^), sma{z), 0) . (C2) 

One can readily check that 

fc"fc„=0, k>^Vf3ka^0, (C3) 

confirming that (CI) is a congruence of null geodesies afhnely parametrized by v. 
The first comoving coordinate 

Z^ = z (C4) 

is trivial: It determines the plane in which the geodesic lies. The second comoving coordinate Z'^ is the coordinate y' 
of the cartesian system {x', y', z) obtained from {x, y, z) by the rotation about the ^-axis by the angle a{z): 

Z"^ = —xsma{z) + y cosa{z) . (C5) 

In the rotated cartesian systems (a;', y', z), the particles move along the x'-axes, with y' = const. The third comoving 
coordinate Z^ is the retarded time U- corresponding to that direction: 

Z^ = U- = t — x' = t — X cos a{z) — y sin a{z) . (C6) 

From Eqs. (C4)-(C6) we obtain the covectors Z'^ '■ 

Z\„ = (0, 0, 0, 1), 

Z"^ = (0, — sina(z), cosq:(z), —xa^z) cosa{z) — ya'{z) sinQ:(z)) , (C7) 
Z^ — (1, ~ cosQ!(2;), — sinQ:(z), xa' (z) sinQ:(z) — ya' (z) cosa(z)) , 

where a' •.= da/dz. It is easy to see that Z'^^a are independent covectors. Since 

ka = { — 1, cosa{z), sina(.2), 0) , (C8) 

the decomposition (5.5) is obtained with the coeflacients 

wi = xa'{z) sin a{z) — y a' {z) cos a{z) , IU2 = , = —1 . (C9) 

One can easily check that iVk are constant along the geodesies, Eq. (5.8). This also follows from Eq. (CI) which 
allows us to write wi in the form wi = a'{zQ) (xq sinQ!(zo) — jyo cosq:(zo)) • Similarly, one can check that k°'Z'',a^O, 
as given by Eq. (4.6). One can also check the fact mentioned in Section 4, that for a twisting congruence all vectors 
Z'^^a are spacelike (except perhaps a set of measure zero). The spacelike character of the vectors Z^^a and Z^_o, is 
evident; for Z^^a we have 

r)°'f^Z^,„Z^^(3 = a'{zf {x sin a{z)-y cos a{z)f. (CIO) 
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The vector Z^ ^ is thus spacehke unless a' = 0. Calculating the twist cu of our congruence (see Eq. (3.24)), we find 

u; = i|a'|. (CU) 

Hence, if the congruence is twisting, all the three vectors Z'^ a are spacelike. When a' = 0, ka = Z^ a , so that the 
congruence is hypersurface orthogonal. 

Instead of the comoving coordinates Z'', one can, of course, use other comoving variables Z*' = Z'^ [Z'^)- Also, one 
can parametrize the geodesies by a label u different from the affine parameter v. When we change the parameterization, 
V = v{u, Z'^), the null vectors are rescaled: 

fc" (7" = ^ = ^k". (C12) 
du ou 

This leads to a new decomposition, namely 

[/„-W2Z2,„+W3Z3,,, (C13) 

where W2 = {dv / du)w2 , W3 = dv/du . As discussed in Section 5, if the congruence is not affinely parametrized, i.e., 
if dv/du 7^ const ^ the coefficients W^ are not necessarily comoving. 

II. The second example will be described only briefiy. It is the familiar ingoing principal null congruence in Kerr 
spacetime. In the ingoing Kerr coordinates {V , r, 9, (p) which generalize the ingoing Eddington-Finkelstein coordinates 
of the Schwarzschild metric, the Kerr metric reads (our notation follows |^): 

ds^ = - (1 - 2Mrp-^)dV^ + 2drdV + p^dO"^ 

+ P'^ [{r^ + a^ f - Aa^ sin^ 6] sin^ Odip^ (C14) 
- 2asin^ 9dp>dr - AaMrp^'^ sin^ 6d(pdV . 

Here, the constant parameters M and a are the mass and angular momentum per unit mass, and the functions A 
and p have the form 

A = - 2Afr + a2, ^ + cos^ 6*. (C15) 

The ingoing null Kerr congruence is given by 

V = const, r = —v, 9 = const, ip = const, (d6) 

where we have absorbed a constant energy parameter into the affine parameter v (cf. |^). The coordinates 
Z-^ := V, Z^ :— 9, Z^ :— (p are clearly comoving. We can easily form the basis vectors Z^ ■ 

Z\^ = (1, 0, 0, 0), 

= (0, 0, 1, 0), (C17) 
Z\c. = (0, 0, 0, 1). 

The covariant components of the tangent null vector fc" — dx" /dv are 

fca = (-1, 0, 0, asin^ei). (C18) 

Its decomposition into the three covectors Z'^ ^ yields the coefficients 

wi — ~1 , ^02=0, W3 — asin^ 9 , (C19) 

which are constant along the null geodesies (C16), in accordance with Eq. (5.8). The covariant metric can be read 
off from Eq. (C14). The norms of the vectors Z'^ are 

g"^Z\aZ\fj = p-^a^ siii^9, 

g''Pz\^Z\p^p-\ (C20) 
g''f'z\^Z\fi^{psin9)-^, 
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where is given by (C15). We see that all the vectors ^ are spacelike as long as a 7^ 0, i.e., when the congruence 
is twisting. The twist w, given by Eq. (3.24), is 

uj = \a COB e\p-'^. (C21) 

The congruence (C16) is twisting even in the flat-space limit of the Kerr metric, obtained by putting M = 0. In 
fact, Eqs. (C20) and (C21) are independent of M . With a = 0, the vector becomes null and the congruence is 
hypersurface orthogonal, ka = —Z^^^ , i-e., nontwisting. 

The comoving coordinates Z^ = 9, Z^ — tp are simple, but they become singular at the axis 6 = and 6 = tt , the 
magnitude of the vector Z^ ^ becoming infinite, Eq. (C20). It is easy, however, to cure this defect by going over to 
another pair of comoving coordinates, Z^ and Z^ , e.g., 

Z^' = sin (fi sin 6, Z^' = cos if sin 9 . (C22) 

Then 

ka = —Z^a + asm9cos(pZ^ — a sin sin (^Z^ , (C23) 

and 

g"'^ Z'^\aZ'^\p = p-^{l - siri^ ip sin"^ 9), 

g°'P Z^\aZ^\l3 = p-^{\ - cos^ sin^ 9) (C24) 

are regular at 6* = and 9 = tt. 
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